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Abstract 

Let 

£ i> • • • j £m be i.i.d. random variables with 

P(Ei = 1) = P(Ei = -1) = 1/2, 

and X m = Y^=i £ i- Let Y m be a normal random variable with the 
same first two moments as that of X m . There is a uniquely deter- 
mined function ^f m such that the distribution of ^> m (Y m ) equals to 
the distribution of X m . Tusnady's inequality states that 

Y 2 

\Vm(Y m )-Y m \< -^ + 1. 

m 

Here we propose a sharpened version of this inequality 
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1 Conjecture 

Let ei, . . . , e m be i.i.d. random variables with 

P(£, = 1) = P(s t = -1) = 1/2, 

and X m = 'Y^ T ° =l £i- Let Y m be a normal random variable with the same 
first two moments as that of X m . Using quantile transformation we can 
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see that there is a uniquely determined function \l/ m such that the dis- 
tribution of ty m (Y m ) equals to the distribution of X m . The central limit 
theorem implies that the function \I/ m is close to the identity for large m. 
A sharp inequality of Tusnady p2] raised certain interest in the literature 

mmmmmmmmmmmm)- 

Let us define the function / on the interval (0, 1) as 

f(x) = y/(l + x) log(l + x) + (l-x) log(l - x), 
set f(0) = 0, f(l) = v/log(4). Let us put 

h — 22 
2 

X k,m — m 
2 

for positive even integers m with k such that m/2 < k < m, and set 

p Km = P{X m >2k-m) = 2- m J2[ m )- 

i=k ^ ' 

Let us define the function Q on the reals as 

Q(x) = — = / e- u2/2 du. 

With those ingredients our conjecture states that 

Q{y/rnf(x k , m )) < p k , m < Q(Vmf(x k -i )m )) 
holds true for y < k < m. Or more sharply 

2(k - 1) - - + 0.8964 < mf-^Q-^pk^/Vm) < 2(k-l) - — + 1.0000 (1) 

holds true with pessimal parameters m — k — 10. It implies that Tusnady's 
inequality is sharpened to 



in 



< 1.1036. 



2 Generalization 

For an arbitrary random variable X let us consider the function on reals 

R(t) = Ee tx 
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restricting ourselves for distributions having finite momentum generators. 
Next we define 

, , \ R'(t) 

m = Wv 

a(x) = t iff ip{t) = x, 

p(x) = R(a(x)) exp(— xa(x)). 

The probability P(^2T =1 Xi > mx) is approximately p(x)~ m if x > EX. 
The function p depends on the distribution of X, it is the Chernoff function 
of X. Let us denote the Chernoff function of the distribution F of X by 
Pf, and the corresponding function for standard normal by pa- The quantile 
transformation between the partial sums of distribution F with Gaussian 
ones resemble us to the equation 

Pf{x) = p G {y) 

having the property that it gives sharp values for any m. Perhaps the error 
term is bounded with a bound depending on the distribution of X. For the 
case symmetrical binomial distribution the error term might be as small as 
that the quantile curve jumps over its limiting function: it is the informal 
explanation of our conjecture. 



3 Numerical Illustration 

The function \l/ m is shown in Figure 1. called "step" for m = 50 with a 
rescaling for random variables 

£ _ -S^-m 1 m 

Sm j Vm ■ 

m m 

The function / is called "limit", for the sequence of step functions goes to 
/ after rescaling. The conjecture comes from the observation that the limit 
function crosses all steps near to their middle. Let us introduce the blow up 
error term 



10 [2k - 1 - m 



^(S(7H)> 



for < k < m/2. In Figure 1. it is labelled as "Delta". With these notations 
([I]) is equivalent with < A*. jm < 1.036. These error terms are shown in 
Figure 2. for 2 < m < 1000. Figure 2. prompts the conjecture that even 
these curves are convergent. We are a bit perplexed: even the inequality 
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< Ai, 2 < 1.036 means that Q(0.723359) < 0.25 < Q(0.6435214). How can 
we prove such an inequality theoretically? 

Acknowledgement: We thank to Peter Harremoes for pointing out a 
mistake in the earlier version of the paper. 
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4 Appendix 



R- program of Figures 1 and 2. 
Q=f unction(p) -qnorm(p) 

G=f unction (x) ((l+x)*log(l+x)+(l-x)*log(l-x))**0.5 
Ginv=f unction(u) { 
GG=f unction(x) G(x)-u 

uniroot (GG , c (0 , 1) , f . lower=-u, f . upper=log(4) " . 5-u , tol=10~-100) 
} 

m=50; k=m/2 

sum=0; divisor=2**m; bin= 
xx=c(l:k+l); yy=c(l:k+l); zz=c(l:k+l); 

for (i in l:k-l){ 

sum=sum+bin 

x=(m-2*i)/m 

y=Q (sum/divisor) / (m** . 5) 

b=Ginv(y)$root 

yy[i+l]=y; xx[i+l]=x 

bin=(m-i)*bin/(i+l) 

zz [i+1] =10* (m-2*i-l-m*b) } 

xx[k+l]=0; yy[k+l]=0; zz[k+l]=0 

kerx=c (0,1.25) ; kery=c (0 , 1 . 15) 

plot(kerx, kery, type="n" ,xlab="eta" , ylab="xi", 

main= " Figure 1 . Quantile transform, its limit and blownup error, m=50") 
for (i in l:k){ 

bb=seq(f rom=yy [i+1] , to=yy[i] , by=0.01) 
cc=bb*0+l; cc=cc*xx[i+l] 

points (bb, cc, type="l" , col="blue", lwd=2)} 
cc=seq(from=0, to=0.999, by=0.001) 
bb=((l+cc)*log(l+cc)+(l-cc)*log(l-cc))**0.5 
points (bb,cc, type="l", col="red", lwd=2) 
points(yy,zz, type="l", col="green", lwd=2) 
legend(locator(l) ,c("Limit" , "Step" , "Delta") , 
lty=c(l,l,l), 

col=c("red" , "blue" , "green")) 
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kerx=c (0,1.25) ; kery=c (0 , 1 . 15) 

plot(kerx, kery, type="n" ,xlab="eta" , ylab="Delta" , 
main="Figure 2. The blownup error") 

for (k in 1 :500){m=2*k; 
sum=0; divisor=2**m; bin=l 
yy=c(l:k+l); zz=c(l:k+l); 
for (i in l:k-l){ 
sum=sum+bin 

y=Q (sum/divisor) / (m** . 5) 

b=Ginv(y)$root 

yy[i+l]=y; 

bin=(m-i)*bin/(i+l) 

zz [i+1] =10* (m-2*i-l-m*b)} 

yy[k+l]=0; zz[k+l]=0 

if (k<100) clr="red" else 

if (k<200) clr="blue" else 

if (k<300) clr="purple" else 

if (k<400) clr="gray" else clr="green" 

points(yy,zz, type="l", col=clr)} 

legend(locator(l) ,c( M 0<m <= 200" , "200<m<=400" , "400<m<=600" , 

"600<m<=800" , "800<m<=1000") , 

lty=c(l,l,l,l,D, 

col=c("red" , "blue" , "purple" , "gray" , "green") ) 
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Figurel. Quantile transform, its limit and blownup error, m=50 
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Figure 2. The blownup error 
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